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Abstract. We consider a class of Jacobi matrices with unbounded coefS- 
cients. This class is known to exhibit a first-order phase transition in the sense 
that, as a parameter is varied, one has purely discrete spectrum below the 
transition point and purely absolutely continuous spectrum above the transi- 
tion point. We determine the spectral type and solution asymptotics at the 
transition point. 
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\0 ■ 1- Introduction 



In this paper we analyze spectral properties of Jacobi matrices, 

/bi ai •••\ 

ai ^2 02 
02 63 as ■ ■ 



(1) J = 

V: : : : -.y 

acting in ^^(Z-(_), and asymptotic properties of solutions to the associated difference 
' equation 

(2) anUn+l + bnUn + an-lUn-l = EUn- 



Recently, there has been interest in the case of unbounded coefficients a„,6„; 
see, for example, PElEllHlElIiniiniII21IIHl- 

Motivated in particular by El El Ej , Janas and Naboko .1^ studied a large 
class of Jacobi matrices with unbounded and periodically modulated entries. To 
be specific, they considered the case a„ = c„/i„, 5„ — dnVn, where {c„} is strictly 
positive and A^-periodic and {d„} is M-periodic. The sequences {/in}, {rn} are 
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unbounded and satisfy a number of conditions. The main example one should have 
in mind is /i„ — Vn — , where < a < 1. Let us consider this special case, that 
is, 

(3) a„ = c„n", 6„ = d„n". 

Since the Carleman condition holds, that is, 

oo 
n=l 

J defines a self-adjoint operator in £^(Z+). The spectral type of J is closely related 
to the location of zero relative to the spectrum of the associated periodic Jacobi 
matrix Jpcr which is given by 

fdi ci •••\ 

Ci ^2 C2 • • 

C2 da C3 • • 



V: : : : 

Clearly, Jpor is if-periodic, where K is the least common multiple of M and N. Its 
characteristic polynomial, d{E), is given by 

K 



d{E) = Tr 



n 







1 



{E - d^)c-^ 

It is well known that the spectrum of Jpor is given by 

f7(Jper) - e M : \d{E)\<2}. 

This set is the union of K bands (non-degenerate closed intervals) whose interiors 
are mutually disjoint. The following was shown in |10j : 

Theorem 1 ( Janas-Naboko) . (a) // |c?(0)| < 2, then the spectrum of J is purely 
absolutely continuous and cr{J) = M. 

(b) //|d(0)| > 2, then the spectrum of J is purely discrete. 

In other words, if zero is not one of the band edges of (t( Jper), one has a complete 
understanding of the spectral type of J. Moreover, Janas and Naboko also described 
the asymptotic behavior of the solutions of compare TO, Theorems 3.1 and 4.2]. 

The question of what happens at transition points, corresponding to d(0) — ±2, 
was left open in UHj. It was pointed out that new methods and ideas would be 
necessary to understand these critical cases. It is our goal here to study this scenario 
in a simple special case. We shall see that even in this simple situation, the analysis 
is already quite involved. Moreover, in this way the main new ideas are more 
transparent. 

We will study the case where the a„ 's and 6„ 's are given by lO and both periods, 
M and N, are small. Specifically, let us consider the case 



M 



5, C?2 



We find 

d{0) = -2 + bb. 

Fix 6 > 0. Then, by Theorem ^ J has purely absolutely continuous spectrum for 
(not too large) 6 > and purely discrete spectrum for 6 < 0. Similarly, for 6 < 0, 
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J has purely absolutely continuous spectrum for 6 < (again, the absolutely value 
should not be too large) and purely discrete spectrum for 6 > 0. Thus, if we fix 
some non-zero value for 6, the case 6 = is the critical case for which Theorem^ 
does not apply. 

It is our goal to study this particular case, that is, we want to determine the 
spectral type and solution asymptotics for 

_ a I, _ / if ^ is odd, 

[0 II n IS even, 

where h ^ Q and < a < 1. Our main results in the case & > (when 6 < 0, one 
has to reflect the energy about zero) are as follows: 

• If 2/3 < a < 1, the spectrum of J is purely absolutely continuous on 
(— cxD,0). Moreover, explicit solution asymptotics are given. These results 
are stated in more detail and proven in Section |21 

• If < a < 1, zero is not an eigenvalue of J. See Section 0] 

• If < a < 1, the spectrum of J is purely discrete in (0, oo). The eigenvalues 
are simple and the n-th eigenvalue, £'„, obeys the bounds 

Ci(6)n" < En < C2{b)n". 

We also provide explicit expressions for the (positive and finite) constants 
Ci (6), 6*2(6). These results and their proofs can be found in Sectional 

This determines the spectral type completely when 2/3 < a < 1. In this case, 
the spectrum looks essentially like 



O'disc(^) 



a'ac(J) 



The condition 2/3<a<lis naturally associated with our method of proof 
(compare, e.g., (|13|l 'l. However, we expect the picture above also when < a < 2/3. 
Thus, we leave the question of proving purely absolutely continuous spectrum on 
the negative energy axis for these values of a as an open problem. Other open 
problems suggested by our work will be discussed in Section |^ 

The spectral analysis in the energy region (—00, 0] is based on an analysis of 
the solutions to the difference equation We obtain asymptotic expressions for 
all solutions corresponding to energies E G (— oo,0]; see Theorems |21 and 0] In 
particular, this determines the asymptotic behavior of the orthogonal polynomials 
associated with the spectral measure of the pair {J, 61) since, by standard theory, 
they solve (0) with a Dirichlet boundary condition, uq = 0. 



2. Preliminaries 
We will study the Jacobi matrix 



(4) 



J = 



ai 




0.2 




02 
&3 as 
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acting in ^^(Z+), where the parameters a„, 6„ are given by 

f> 7 f bn" if n is odd, 

(5) an=n , bn = < „ . 

[ II n IS even. 

Here, b ^ 0. Note that if J{a„}.{b„} denotes the Jacobi matrix corresponding to the 
sequences {a„} and {bn}, and U denotes the unitary transformation of given by 
(C/V')„ = i-ir^n, then 

U J{a„},{b„}U = -J{a„},{-b„}- 

We will therefore restrict our attention in what follows to the case 6 > 0. Consider 
solutions of the difference equation 

(6) anUn+i + b„Un + a„_iU„_i = EUn- 

Defining 

Un-l 



Un = 

the recursion © is equivalent to 

Un+l — TnUrn 



where 



Let 



T = 



1 

a„-i E-b„ 



Bn — T2nT2n-l, 

so that 

(7) (72„+i = (S„ X ... X 



We have 

Bn — 








1 


-(1 


1 \a 
2n) 


E 
(2n)" 


-1 




1 







(2n)' 



1 

_ _J_\a E-b{2n-ir 

^ 2n-l> (2n-l)° 

Q E ' 



{2nY \ -E -bE 














a 1 



0{n 



We see from {TJ that we should study the product of the form i?„ x ... x Bi. 
However, we shall study a slightly different sequence of products whose usefulness 
will become clearer below. To define this auxiliary problem, we have to introduce 
a few matrices. Let 



Cn — 

and 



I -b\ If bE + E/{2b) \ 1/00 



1 J (2n)" V ~E/2 J ' 2n\ ~a ab 

. 1 \ 1/00 

Bn — 



-I 2 J (2n)" \ bE J n \ a -a 



Lemma 2.1. We have 

CnBnCn' = -B„ + 0(n~2"). 

Proof. This is tedious but straightforward. □ 
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3. The Absolutely Continuous Spectrum and Solution Asymptotics 

AT Negative Energies 

Motivated by Q and Lemma |2. II we will study asymptotics for the auxiliary 
problem 

(8) K+i = (S„ X ••• X 

This problem is a more general version of the one studied in jOj Section 3]. There, 
products of matrices of the form 

, 1 \ , 1 / 

tin 



-1 2 ) n\l A 

were studied. We will employ a similar strategy to find solution asymtotics for the 
more general problem at hand. This explains why we introduced the matrices Bn 
and derived Lemma |2. II in the previous section. 

Theorem 2. Suppose that 2/3 < a < 1 and b > 0. Then, for every E < Q, ^ has 
two linearly independent solutions 

(9) - 
with asymptotics given by 



(10) vt = n-"/4 exp ( ±*^-^ • ) (1 + o(l)) asn^^. 



.V^bE 

2"/2 

Proof. We make the ansatz 

Zn = n'' exp (^An^^ 

and define the matrix 

g / Zn — 1 Zn—1 

\ Zn Z^n 

Our goal is to choose 7, A, S such that 

(11) S-l^BnSn = I + Rn, {\\Rn\\} ^ 

Consequently, an arbitrary non-trivial solution of (jSJ has the form Vn = SnWn, 
where Wn is a sequence of vectors which tends to a non-zero vector W. Since 111|) 
will be shown to hold if we let 



(12) ^--4'^=W^(T3|)'^ = l-2' 

the assertion of the theorem then follows immediately. 

Let us consider the matrix S^^iBnSn- It is readily checked that 



where 



and 
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where 

Putting this together, we obtain 



SnliBnSn — (detS'„+i) ^ 



b„ 

-aZ 



where 



an = ZnZn-i + z.^z^+i - 2|z„p + [-hE\zn\^] + i [-az„z„_i + a|z„| 

K = ZnZn-l + ZnZn+l - + [-^Ez"^ + i [-Q!2;„Z„_i + aZ^] . 

We consider first the ofF-diagonal elements of S^^iBnSn- 



{detSn+iy^bn 



-2- 



(2n)° 



n L 



For the denominator in the last expression, we find 

Zn+iz^z'"^ - z;iqrrz^^ = 2A5n^^^{l +o(l)), 

provided we choose 5 as in H12|l (which implies 1/2 < 6 < 2/3). 

For the numerator, we find after lengthy but straightforward calculations. 



bE 
(2n)° 



5-2 



where 



6£; 



25-2 



X = A5{5 - 1) + 2-1 AS + F = {ASy - — 



If 7, A, (5 are chosen as in (|12|) . then X = y = and hence 
(13) (det5„+i)-i6„ = 0(n3^-3) ^ ^ 

which is summable by our assumption on a. 

Next, we consider the diagonal elements of S~]^^BnSn'- 

|a„ — detS'„_|_i| = z„z„_i + z„+iz„ — 2|z„p — ^^^^j'j — h ^ [— z„z„_i + |z„|^ 

^n^n — 1 ^ ^7i+l^n ^jz^^j ~(2n)^^ "n [ ^n^n — l ~t~ |^n| 
— ^n^n — 1 ^ ^-^n ^2n)" n" [ ^n^n — 1 "1" ^7^] 

and hence | (det 5„+i)^^a„ ~ l| = |(det S'„+i)^^6„|. This shows that ((TT|l holds, 
concluding the proof. □ 

The special case a — 1 deserves an additional remark. Theorem |5] above gives 
the asymptotics 



(14) 



"1/" exp [±W~2hE'n^ (1 + o(l)) as n oo 
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for a pair of linearly independent solutions of the auxiliary problem (jHJ . This can 
also be shown by an application of the Birkhoff-Adams Theorem; see Theo- 
rem 8.36]. The latter result concerns second-order difference equations of the form 



(15) x{n + 2)+pi{n)x{n + 1) + p2{n)x{n) = 0, 

where pi{n) andp2('^) have asymptotic expansions 

OG OO , 



with do ^ 0. In our situation, we have 



bE 

Co = -2, ci = 1 do = 1, di 



= -1. 



We have to compute the roots Ai, A2 of the characteristic equation A^ + coA + do = 
and find Ai = A2 = 1. Moreover, if i? 7^ 0, we have 2di ^ cqCi. Thus, part (b) of 
O Theorem 8.36] tells us that there are two linearly independent solutions xi{n), 
X2{n) of H15|l whose asymptotics are given by H14|l . (Note, however, that there is a 
misprint in formula (8.6.7) of "F: rf' should be replaced by n".) 

Theorem 3. Suppose that 2/3 < a < 1 and b > 0. Then, for every E < 0, ^ has 
two linearly independent solutions 



U. 



2n+l 



^2n+l 



with asymptotics given by 



■0(1)), 



where 

T=l( ' 
b\ -I I 

and are given by (O and (|10|) . Moreover, the spectrum of J is purely absolutely 
continuous on (— cxd,0). 

Proof. By standard results of asymptotic analysis (centered around Levinson's The- 
orem; compare Section 8.3] and [S]), Lemma [2. II implies that the solutions of 
JTI) behave asymptotically like the solutions of the system 



U2n+1 — 



( — C„ BnCn 



X (-CfiBiCi) 



Since CmC,„_;^ 



/ + 0(m ^), we may as well consider the system 



U2n+1 = (-1)"C„ 



Now the first assertion is an immediate consequence of Theorem 12 

To prove that the spectrum of J is purely absolutely continuous on (—00, 0), it 
suffices to show that for every E < 0, there does not exist a subordinate solution of 
(jSJ in the sense of Gilbert-Pearson (cf. 0; see also ^1^]). Thus, we claim that 
for every pair of non-trivial solutions u^^\ u^"^^ of lO, we have 

2- 



(16) 



lim sup 

N^oo 



N 



,(1) 



N 



,(2) 



> 0. 
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Since every solution u of © is given by a linear combination of the two solutions 
whose asymptotic behavior is given by (-l)"rv;±(l+o(l)), there are (u-dependent) 
positive constants Ci,C2 such that 

N 
ra=l 

From this, Hlt)|) follows immediately. □ 



4. Solution Asymptotics at Zero Energy 



In this section we study the asymptotics of the matrix product i?„ x • • • x Bi, 
and hence the asymptotics of solutions, for E = 0. We show in particular that 
there are no non-trivial solutions in £'^. 

Theorem 4. For zero energy, E = 0, the difference equation (j^J has two linearly 
independent solutions Un\ j — 1,2, with asymptotics given by 



,(1) 



(-l)"n-t(l + o(l)) 




^2n+l 



(-1)" [feni-T + 0(n-t logn)] (1 + o(l)) 
(1 + 0(1)) 



In particular, zero is not an eigenvalue of J. 
Proof. We have 



Bn — 



1 \ / 1 

-(1-2^)" ; 1, -{i-^,r -b 

-1 -b 

-1 



— 

2n " 



with bounded sequences {r„}, {s„}. 
That is, 



5„ = -Mfc + — / 
2n 



\ , fib 
0(n-2) ) ' ^here M, = I ^ ^ 



Therefore, 

i?„ X • • • X i3i = (-1)"M," (Af-"(-i3„)A4"-i) X • • • X {M-\-B,)M^) 
Since 



M^^(-B,)Af;?-i=/-- 



2j V 1 



o(r') '' 



we obtain 



B„ X ••• X Bi = (-1)' 



1 bn 
1 



n 
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where 



1 + O0-2) 

1 + 0(i-2) 



o(r^) 





Ci+o{l) O(logn) 
C2 + 0(1) 

Here, Ci,C2 are suitable non-zero constants. 
On the other hand, 



with a suitable constant C3 ^ 0. 

Putting everything together, wc find that 

B X xB - r n"r,-f 1^ ^1^3+0(1) C2C36n + 0(logn) 
The assertion follows. 

5. The Discrete Spectrum 



□ 



In this section we study positive energies and prove that the spectrum in (0, 00) 
is purely discrete and accumulates only at 00. Moreover, we provide upper and 
lower bounds for eigenvalues. 

Fix the parameters 6 > and < a < 1. For a > 0, we define 



span 



|<52„-i:l<n<-(^- 



1/q ~ 



Lemma 5.1. For a > and V' e H^a^ n D{J), we have \\{J - a/)V'|| > a\\ 
Proof. Consider the even and odd subspaces of 

= {/ e /2„-i = Vn e Z+}, il = {fef: /2„ = OVn e z+} 
and the respective orthogonal projections, P^ and Pq- We write a vector ip G i"^ as 

V' = V'e + V'o = Pc^ + PoV'- 

Let S be the shift, given by {S^)n = ^n+i, and let A and B be diagonal matrices 
with An^n = dn and -B„,n = 6„, respectively. Then we can write the Jacobi matrix 
J in the form 

If we write 



J^SA + AS* 



B. 



Po{J- 


aI)Po 


Po{J- 


aI)Po 


Pe{J- 


aI)Po 


Pe{J- 


aI)Pe 



and use 



we obtain 



J-aI = 



BP^ = 0, Po{SA + AS*)Po = 0, P^{SA + A5*)Pe = 0, 



(P - ar)P, 


PoJaP. 


PeJAPo 


-aP^ 



J-al 

Here, Ja denotes the matrix SA + AS* . This yields 
{J-alf-a'^I 



(P2 _ 2aP)Po + (Pe JaPo)* (PeJAPo) 


{B - aI)PoJAPe - aPoJAPe 


P^JaPo{B - aI)Po - uP^JaPo 


{PoJAPeTiPoJAPe) 
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We want to show that 

(17) ([(J - a/)2 - a^I]ij,i^) > for e T^i^). 
Observe first that 

(18) (^2 - 2aS)Po > on Hi^), 
since 

(6(2n ^ l)"f - 2ab{2n - 1)" > ^ 2n - 1 > 
Now, let tp € Tia n D{J) be given and write 

Ipe = Pci>, V'o = Polp, (t>c = Pc JaV'o, (/-o = PoJa^c- 

Then 

([( J- a/)2 - a2/]V^, V) = ((B' - 2aB)^o, i'o) + ||0c||' + ||(/)o|P + 2Re((B- 2a)0o, ^o). 
In view of (|18|l . (|17|l follows from this once we show that 

(19) |((i?-2a)0o>o)| < ^ [((B2-2ai?)Vo,V'o) + ||0c|P + ||'^o|P] ■ 
Since 

\{iB - 2a)0o>o)| < ||0o|| • ||(S - 2a)^o|| < ^||'/'o||' + ^ ■ \\{B - 2a)V^|P, 
((T^ follows from 

(20) i-||(B-2a)V'o|P< ((S2-2aB)^o,V^o). 
But (pn|) is a consequence of 

0< ((^2 _4a2)^„,^„) 

which in turn follows immediately from the assumption ip G Tla ■ □ 

The following theorem describes the positive spectrum of J, which turns out 
to be discrete in [0, cxd). It is therefore interesting to study the distribution of 
eigenvalues 

< Ei{b) < E2{b) < ■■■ 

Note that all eigenvalues must be simple. Let N{E) denote the number of eigen- 
values of J in the interval (0, E). 

Theorem 5. Let < a < 1 and b> 0. 

(a) o'(J) n (0,oo) is discrete and can accumulate only at oo, not at 0. 

(b) The eigenvalue counting function N obeys the estimate 

«(.)a(f)*- 

(c) The n-th eigenvalue, En{b), obeys the lower bound 

En{b) > 2"6n". 

(d) Ifb>V6, then 

22+a 

EJb) < -— —bn". 
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(e) Ifh<V&, then 

for some constant C > 0. 

Proof, (a) and (b) are immediate consequences of Lemma |5.1I On 7ii^\ we have 
(J — a/)^ > a^, and hence the dimension of the range of the spectral projection of 
J onto the interval (0, 2a) is bounded by 

(c) follows from (b): For £ > 0, we have 

1 /i?„(&) + e^^/" 



n = NiEn{b)+s) < 2 
and hence 

En{b)+e> 2"6n". 

Now let £ go to zero. This yields the claimed lower bound for En{b). 

In order to prove the upper bounds in (d) and (e), we note the following. If, for 
n > 1 and a > 0, we can find test functions /*^^\ . . . , /'■"•' with disjoint supports 
obeying the estimate 

(21) |l(J-a/)/^^^||' <a'||,/^''^||' 

for 1 < fc < n, then E„{b) < 2a. 

We first prove the estimate in (d). Our test functions will be of the form S2m-i, 
where m is chosen such that 

|6(2m-ir-a|<^. 

Equivalently, 

For these test functions, we have 

||(J - aI)S2,n-if - a'||52m-i||' < 2(2m - if" + \b{2m - 1)" - af - a" 

<26^ + T-^- 

If 6 > -s/G, this shows that ||(J — al)52m-i\\'^ — a^||<52m-i|P is negative. Note that 
the number of test functions we obtain this way, n, is restricted by the condition 



1 

n< - 
- 2 



26 



f-V 

\2bJ 



Thus, if we choose 

2n 

a> 2b 



^31/" - I 

we can find n test functions with disjoint supports obeying H21|l. This yields the 
estimate in (d). 



12 



D. DAMANIK, S. NABOKO 



Let US turn to part (e). Again, our test functions will be supported on odd 
sites 2m — 1, subject to H22|l . We partition the interval in (|22|l into n subintervals 
Ji , . . . , J„ of equal size, roughly given by 



A. 



2b) \2b) 



Since we need this to be at least 2, we get a first condition for a: 

(23) A„ > 2. 

We will return to this condition below. 

The function /^'^^ will be supported on the odd sites within Jk and on these 
sites, it alternates between the values 1 and —1. Our goal is to prove the estimate 
By definition of f'-''\ we have that 



A. 



Using the fact that Jj. lies within the interval H22I) . we obtain 



|(J-a/)/Wf < A„. sup |6(2m-l)^ 

2m-le,/fc 



const • sup (2m) " + 

2m-leJk 



J2 |(2m+ 1)" - (2m- 1)"| 



2m-leJk 

< An h const ' 

4 



1/a 



const -An sup (2m-l)2(°-i) 

2m-leJfc 



with suitable uniform constants. Thus, in order to satisfy H21|) . we need 



(24) const • 

Let us first prove 
(25) 



3a . 

— + A„ sup (2m - 1 



2b 



,2(q-1) 



2m-leJfc 



< A„-. 



const ■ 



This inequality is equivalent to 
(26) ^ 



'3^1/" 



'1^1/a 



< A. 



^l-2a^Q < a. 



Note that the factor in front of b^ ^" is uniformly bounded from above by some 
C > for all a € (0,1]. 

Since, by 2m - 1 > (^)^^" on Jk, we have 

2(q-1) 



A„ sup (2m- l)2("-i' < A^ 

2m-leJfc 

Thus, to complement (|25|l . it suffices to show 



1/a 



< A, 



(27) const • A^ 
This inequality is equivalent to 

(28) (8-const)"/2(25)i-" < a. 
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Since we already have the condition a> C ■ fi-om H26|l . (|28|1 is automatically 

satisfied if C obeys, in addition to the property above, 

C> max 8-const"/22i-"(^/6)". 

0<a<l 

(Here we used b < \/6.) 

Finally, we need to satisfy (|23|l . Using a > Cb^~'^"n°' , one checks that A„ > 2 
holds as long as 



C 



(1) 



V._(,)l/c 



> 1. 



2" (76)2" 

(Here we used again that b<y/6.) This can again be satisfied uniformly in a e (0, 1] 
by choosing C large enough. 

Putting everything together, if C is so large that it satisfies the three properties 
found above, then for 

a > Cb^^^°'n°', 

we can find n test functions f^^\ . . . , /'"^ obeying H21|l . □ 

6. Open Problems 

In this section we list a number of questions and directions for future research 
that are suggested by our and previous work. 

(1) More general entries: We expect that one can relax the assumption = 
r„ = n°' and perform a similar analysis, akin to 10 , for more general 
unbounded sequences. 

(2) The n-dependence of i?„(6): In fact, a stronger statement on the behavior 
of En{b) for large n would be desirable. Is it possible to find a function 
F{b) such that En{b) ^ F{b)n°' (in the sense that En{b)F{b)-^n-°' 1 as 
n oo)? 

(3) The 6-dependence of E„{b): Some improved estimates on the dependence 
of En{b) on b would be of interest; especially in the case a S (1/2, 1) for 
which the statement in Theorem 5.(e) does not appear to be optimal. 

(4) Eigenfunction asymptotics for positive energies: We established explicit 
asymptotics of the solutions to the associated difference equation for all 
energies in (— oo,0]. Is it possible to prove similar results for energies in 
(0, oo)? We would expect the formula from Theorem |21 to hold also for 
positive energies. 

(5) More detailed questions as b crosses zero: Combining the results of JHl and 
the present paper, we obtain a rather detailed picture of the spectrum and 
the spectral type as b makes a transition through the critical value 6 = 
(cf. the discussion at the end of Section 1). This suggests that even more 
detailed questions could be addressed. For example, what happens to the 
spectral function as b passes through zero? Do we see some concentration 
at the eventual eigenvalues? 

References 

[1] J. Dombrowski, Absolutely continuous measures for systems of orthogonal polynomials with 
unbounded recurrence coefficients, Constr. Approx. 8 (1992), 161-167 

[2] J. Dombrowski and S. Pedersen, Absolute continuity for unbounded Jacobi matrices with 
constant row sums, J. Math. Anal. Appl. 267 (2002), 695-713 



14 D. DAMANIK, S. NABOKO 

[3] J. Dombrowski and S. Pedcrscn, Spectral transition parameters for a class of Jacobi matrices, 
Studia Math. 152 (2002), 217-229 

[4] J. Edward, Spectra of Jacobi matrices, differential equations on the circle, and the su{l, 1) 
Lie Algebra, SIAM J. Math. Anal. 34 (1993), 824-831 

[5] S. N. Elaydi, An Introduction to Difference Equations, 2nd edition, Undergraduate Texts in 
Mathematics, Springer- Verlag, New York (1999) 

[6] M. S. P. Eastham, The Asymptotic Solution of Linear Differential Systems. Applications of 
the Levinson Theorem, London Mathematical Society Monographs. New Series, 4, Oxford 
Science Publications, The Clarendon Press, Oxford University Press, New York (1989) 

[7] D. J. Gilbert and D. B. Pearson, On subordinacy and analysis of the spectrum of one- 
dimensional Sclirodinger operators, J. Math. Anal. Appl. 128 (1987), 30-56 

[8] J. Janas and M. Moszyriski, Spectral properties of Jacobi matrices by asymptotic analysis, 
J. Approx. Theory 120 (2003), 309-336 

[9] J. Janas and S. Naboko, Spectral properties of selfadjoint Jacobi matrices coming from birth 
and death processes, in Recent advances in operator theory and related topics {Szeged, 1999), 
Oper. Theory Adv. Appl. 127, Birkhauser, Basel (2001), 387-397 
[10] J. Janas and S. Naboko, Spectral analysis of selfadjoint Jacobi matrices with periodically 

modulated entries, J. Fund. Anal. 191 (2002), 318-342 
[11] ,1. Janas and S. Naboko, Criteria for semiboundedness in a class of unbounded Jacobi opera- 
tors (Russian), Algebra i Analiz 14 (2002), 158—168; ; translation in St. Petersburg Math. J. 
14 (2003), 479-485 

[12] J. Janas, S. Naboko, and G. Stolz, Spectral theory for a class of periodically perturbed 
unbounded Jacobi matrices: elementary methods, J. Comput. Appl. Matli. 171 (2004), 265— 
276 

[13] S. Jitomirskaya and Y. Last, Power-law subordinacy ajid singular spectra. L Half-line opera- 
tors, Acta Math. 183 (1999), 171-189 

[14] S. Khan and D. B. Pearson, Subordinacy and spectral theory for infinite matrices, Helv. 
Phys. Acta 65 (1992), 505-527 

[15] C. F. Lo, K. L. Liu, and K. M. Ng, The multiquantum intensity-dependent Jaynes-Cummings 
model with the counterrotating terms, Physica A 265 (1999), 557-564 

[16] C. F. Lo. K. L. Liu, and K. M. Ng, Exact cigenstatcs of the intensity-dependent Jaynes- 
Cummings model with the counter-rotating term, Physica A 275 (2000), 463—474 

[17] D. R. Masson and J. Repka, Spectral theory of Jacobi matrices in ^^(Z) and the su(l, 1) Lie 
algebra, SIAM J. Math. Anal. 22 (1991), 1131-1146 

[18] M. Moszynski, Spectral properties of some Jacobi matrices with double weights, J. Math. 
Anal. Appl. 280 (2003), 400-412 



